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PARTICLE  MOTIONS  IN  A  MAGNETIC  FIELD 

Introduction 

The  motion  of  a  charged  particle  in  the  earth’s  magnetic  field 
has  long  been  of  interest  to  mathematicians  and  physicists  in  con¬ 
nection  with  the  study  of  the  pe'.ar  aurora  and  cosmic  rays .  The 

n 

mathematical  formulation  of  this  problem  was  given. by  Stormer  as 
early  as  1907J  it  is  often  referred  to  as  Stormer’ s  problem. 
Recently,  this  problem  received  renewed  significance  with  the  dis¬ 
covery  of  the  Van  Allen  radiation  belt,  (see  Dragt  [6]).  This  is 
a  region  in  space  that  consists  of  electrically  charged  particles, 
which  are  assumed  to  be  trapped  by  the  earth's  magnetic  field. 

Some  of  these  particles  were  observed  to  have  a  lifetime  of  several 
years.  The  purpose  of  this  paper  is  to  rigorously  establish  the 
theory  of  almost  periodic  motions  for  the  Stormer  problem,  exhibit¬ 
ing  thereby  the  trapping  of  charged  particles  as  observed  in  the 
Van  Allen  belt.  An  additional  feature  of  the  theory  we  shall  develop 
is  that  it  can  easily  be  generalized  to  any  rotationally  symmetric 
"mirror  field". 

The  trajectory  of  a  particle  in  a  magnetic  field  is  generally 
very  complicated  and  must  be  obtained  by  numerical  integration  of 
the  differential  equations  of  motion.  In  uhe  special  case  of  a 
uniform  static  magnetic  field  B,  the  trajectories  can  be  obtained 
explicitly.  -As  is  well  known,  ill  particles  gyrate  in  a  helix 
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about  the  magnetic  field  lines  (see  Figure  1). 


Figure  i.  Particle  motions  in  a  constant  magnetic  field. 

If  m  denotes  the  mass  of  the  particle,  q  its  charge,  v^  the 
velocity  perpendicular  to  the  magnetic  field,  and  B  the  magnitude 
of  the  magnetic  field,  then  the  quantities, 

rav^ 

M  =  W  * 
rav 

a=  qS' 

are  constant  along  any  orbit.  M  is  called  the  magnetic  moment  of 
the  particle,  and  "a"  its  radius  of  gyration. 

Many  mathematicians  have  concerned  themselves  with  the  motion 
of  a  charged  particle  in  a  slowly  varying  magnetic  field.  A  slowly 
varying  magnetic  field  is  a  field  which  varies  slowly  in  space  and 
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time  -  that  is,  slowly  compared  with  the  gyration  radius  and  period. 
Essentially  this  means  that  in  the  course  ..of  one  gyration  about  a 
magnetic  field  line,  the  particle  sees  an  approximately  constant 
field.  In  a  slowly  varying  field  the  particle  moves  approximately 
in  a  circle  whose  center  drifts  slowly  across  the  lines  of  force 
mid  moves  rapidly  along  the  lines.  This  is  the  so-called  "guiding 
center"  or  "adiabatic"  approximation.  It  was  shown  by  Alfven  [1] 
that  the  magnetic  moment  is  an  adiabatic  invariant  in  a  slowly  vary¬ 
ing  l'ioldj  that  is  to  say,  it  is  constant  to  first  order  in  the 
radius  of  gyration.  This  result  is  of  extreme  importance  in  plasma 
physics,  where  one  is  interested  in  confining  charged  particles  in 
a  bounded  region.  Suppose,  for  example,  that  the  magnetic  field  is 
a  convex  function  along  the  lines  of  force.  A  particle  moving  a- 
J.on.g  a  Line  of  force  will  be  "reflected"  backwards  at  the  point  Pq 
defined  by 


M3( P  )  =  E 
c 


whore  E  is  the  energy  of  the  particle.  Thus,  to  first  order,  the 
guiding  center  of  a  particle  oscillates  periodically  along  a  line 
of  force,  between  two  "mirror"  point;:.  In  this  case  it  has  been 
shown  (Northrop  [12])  that  the  quantity 


J  »  /  P..ds 
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is  also  an  adiabatic  invariant,  where  P„  is  the  guiding  center  j 

momentum  parallel  to  the  lines  of  force,  and  the  integral  is  taken 

over  a  complete  oscillation  from  one  mirror  point  to  the  other  and  J 

back  again.  J  is  usually  referred  to  as  the  longitudinal  adiabatic  j 

invariant. 

However,  for  virtually  every  prospective  device  for  the  pro-  ! 

duction  of  useful  energy  from  controlled  thermonuclear  fission,  it 
was  seen  that  the  requirement  that  the  particle  remain  confined  for 
periods  of  time  encompassing  many  millions  of  gyrations  could  gen¬ 
erally  be  met  only  if  the  magnetic  moment  were  constant  to  a  much 
higher  order.  In  1955,  Hellwig  [9]  proved  the  constancy  of  the 
magnetic  moment  to  second  order  in  the  radius  of  gyration,  and  in 
1957,  Kruskal  [10J  proved  the  constancy  to  all  orders.  Finally, 

Gardner  [  7]  showed  the  constancy  of  the  longitudinal  adiabatic  in¬ 
variant  to  all  orders.  Moreover,  Gardner  presented  a  general  method 
to  obtain  formal  asymptotic  expansions  for  all  the  adiabatic  in¬ 
variants.  The  main  idea  of  this  paper  is  to  show  that  the  phase 
space  of  a  particle  moving  under  the  influence  of  the  earth' s 
magnetic  field  contains  a  region  where  series,  ialogous  to  the 
formal  expansions  of  Gardner,  are  actually  convergent  expansions. 

This  will  be  accomplished  by  using  a  theorem  of  J,  Moser  which 
guarantees  the  existence  of  almost  periodic  solutions  of  the  differ¬ 
ential  equations  of  motion.*  In  this  manner  we  will  show  that 

Gardner  [8],  in  1>*>2,  announced  a  result  like  ours  for  particle  l 

trajectories  in  a  "mirror"  field.  To  the  author's  knowledge,  a 

proof  of  this  result  was  never  published  by  Gardner.  j 

I 
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particles  which  are  adiabatic  silly  trapped  are,  in  fact,  rigorously 
trapped  for  all  time.  This  possibility  was  first  pointed  out  by 
Arnold  [2], 

The  author  wishes  to  express  his  deepest  gratitude  to  his  thesis 
advisor,  Professor  Jurgen  K.  Moser,  for  his  many  helpful  hints  and 
suggestions  and  above  all  for  his  patience  and  understanding  while 
this  papti  was  being  written. 

1.  The  Stormer  Problem 

The  earth' s  magnetic  field  is  assumed  here  to  be  equivalent  to 
the  field  produced  by  a  magnetic  dipole  situated  at  the  center  of 
the  earth.  Such  a  field  can  be  described  in  cylindrical  coordinates 
p,z,$  by  the  equations 

B*  =  curl  A* 

d.i) 

r- 

B  =  |  B)  =  £L(l+3  sin2X)l/2 
r' 

(see  Figure  2),  where  M  is  the  moment  of  the  magnetic  dipole, 
which  points  in  the  negative  z  direction,  and  £  is  a  unit 
vector  in  the  *  direction.  The  plane  X  »  C  is  the  equatorial 
plane,  and  the  magnetic  lines  of  force  are  given  by 


*1 
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Figure  2. 


From  the  previous  discussion  one  would  intuitively  expect  that 
the  particles  with  small  energy  will  gyrate  about  the  guiding  field 
line  with  the  (so-called  cyclotron)  frequency 


u>  * 


qB 

m** 


» 


whore  vi  and  m  denote  the  charge  and  moss  of  the  particle.  More¬ 
over,  since  the  field  3  is  a  convex  function  along  a  line  of  force, 
we  vc  .id  expect  that  the  particic,  as  it  moves  into  regions  of 
stronger  field  at  higher  latitudes,  will  be  reflected  back  toward 
the  equator  by  converging  lines  of  force.  To  what  extent  this  is 
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true  will  be  discussed  in  the  following  sections. 

To  write  the  differential  equations  of  motion  for  the  Stormer 
problem,  it  is  most  convenient  to  employ  a  canonical  formulation 
described  by  the  Hamiltonian 


(1.3) 


where 


Pp  = 

Pz  =  mi 

2  • 

p,  -  mp  9  +  qpA. 

9 

Since  H  is  independent  of  time,  the  energy 

H  ioV*"  a  E 

is  a  constant  of  the  motion.  A  second  integral  of  the  motion  is 
obtained  by  noting  that  H  is  Independent  of  the  angle  9.  Hence 
the  canonical  angular  momentum 


p0  -  iMr,  (i.i) 

where  F  is  defined  by  this  equation,  is  a  constant  of  the  motion. 


(Xfte  integration  constant  T  had  the  dimensions  of  a  reciprocal 
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length.)  The  three  dimensional  problem  is  now  reduced  to  the 
simpler  problem  of  finding  the  two-dimensional  motion  of  a  particle 
in  the  p  -  z  plane  under  the  influence  of  the  potential 

*<’.■>  =  4r<5  -  h?  ■  cw) 

r 

Once  p(t)  and  2(t)  have  been  found,  </>(t)  is  then  determined  by 
integrating  the  equation 


which  yields 

t 

<t>{ t)  =  m  +  i  f  (Sff  -  .  (1.6) 

The  sign  of  T  plays  a  crucial  role  in  determining  the  general 
properties  of  trajectories.  The  radial  derivative  of  V  is  given. 

by 


r.W  =  Z±.l£(l 
m  P 


P 

3 


2P. 

~y 


(1.7) 


which  is  strictly  less  than  zero  for  T  negativf  .  A  negative 
radial  derivative  for  the  potential  corresponds  to  a  repulsive 
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radial  force,  since  -r«  V  is  the  component  of  the  force  in  the 
radial  direction.  Hence  all  trajectories  characterized  by  a 
negative  r  must  extend  to  infinity  and  cannot  be  trapped.  In 
addition  the  particle  is  restricted  to  lie  in  the  region  V(p,z) 
E.  This  region  is  indicated  in  Figure  3. 


Figure  3.  The  region  V  §  E  for  T  <  0 


VII 


Note  also  that  no  orbits  extend  into  the  dipole  (r  =  0)  for  T 
negative. 

The  situation  is  very  similar  when  P  =  0.  For  p  unequal 
to  zero  the  radial  derivative  of  V  is  again  negative.  However, 
p  =  0  is  a  solution  of  the  equations  of  motion.  Hence  the  tra¬ 
jectory 

z(t)  =  V2mE  t  +  z  ,  z  <0 

runs  into  the  dipole  from  below  the  equator,  and 

z(t)  =  - V§nE  t  +  z  ,  Zq  >  0 

is  a  trajectory  running  into  the  dipole  from  above  the  equator. 
All  orbits  starting  in  the  shaded  region  of  Figure  4  must  extend 
to  infinity. 


Figure  4.  The  region  V  IE  for  T  =  0 


For  the  study  of  bounded  trajectories,  therefore,  we  restrict 
ourselves  to  the  case  r  >  0.  It  is  convenient  at  this  point  to 
introduce  the  dimensionless  variables 


(1.8) 
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It  is  easily  seen  that  the  equations  of  motion  for  these  dimension¬ 
less  variables  are  derived  from  the  new  Hamiltonian 


(1.9) 


where  we  have  omitted  the  primes  for  convenience.  In  this  system 
of  units  the  particle  has  the  dimensionless  velocity 


W 


o 


i 


where 


4  1  qM.2  4 

7i  a  TSW  r 


(1.9)' 


The  dimensionless  constant  f  is  that  used  by  Stormer  [14], 
Note  that  the  angular  momentum  r  is  now  normalized  to  one. 
The  potential 


V(P,«)  a  -  -r)2  (1.10) 

r 

vanishes  along  the  curve  r  =  cos  ^  and  is  positive  elsewhere, 

2 

(The  line  of  force  r  =  cos  X  corresponds  in  our  old  coordinates 

1  2 

to  the  line  of  force  r  «  T  cos  X. )  Since  the  Hamiltonian  H  of 
(1.9)  is  a  constant  of  the  motion,  the  particle  is  restricted  to 


♦ 
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lie  in  the  region  0  §  V  §  H.  This  region  assumes  three  different 
forms  depending  on  whether  H  is  less  than,  equal  to,  or  greater 
than  1/32. 


Figure  5. 


Allowed  region  V(p,z)  S  H 


for 


H  < 


1 

52* 


Figure  6.  Allowed  region  V(P,e)  >  1/32. 
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Figure  7.  Allowed  region  V(P,x)  ^  H  for  H  =  . 

Fran  Figure  5  we  see  that  any  trajectory  starting  in  the  oval 

2 

like  region  surrounding  the  curve  V  =  0  (r  =  cos  X),  with  initial 

energy  less  than  1/32  can  never  leave  this  region  for  otherwise 

it  would  encounter  larger  values  of  V,  The  almost  periodic 

motions  we  shall  find  will  all  lie  in  this  oval  like  region,  where 

the  value  of  H  will  he  very  small.  These  solutions  will  gyrate 

o 

about  the  line  of  force  r  *  cos'Ti  and  oscillate  back  and  forth 
across  the  equator.  Furthermore,  we  shall  show  that  these  motions 
can  penetrate  arbitrarily  close  to  the  dipole,  a  result  which  was 
somewhat  unexpected. 

One  cannot  expect  to  trap  particles  with  H  >  l^S,  since  the 
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region  V  §  H  extends  continuously  to  infinity.  However,  for  just 
this  reason  these  solutions  are  important.  Namely,  a  trajectory 
cannot  extend  into  the  dipole  from  infinity  miles s  H  >  1/52.  Such 
trajectories  play  a  role  in  the  theory  of  the  polar  aurora  (Stormer 

[14]). 

Unfortunately,  there  axe  no  further  known  constants  of  the 
motion,  so  that  the  system  of  equations  derived  from  the  Hamiltonian 
(1,9)  is  as  simple  a  system  as  one  can  achieve.  In  general,  it  has 
no  known  explicit  solutions.  The  equations  can,  however,  be  solved 
in  terms  of  elliptic  functions  for  the  special  initial  conditions 

z  =  z  =  0 

in  which  case  the  orbit  is  confined  to  the  equatorial  plane  (since 
the  magnetic  field  is  perpendicular  to  the  equatorial  plane,  and 
the  force  qv  xf  is  perpendicular  to  B).  The  general  properties 
of  all  equatorial  orbits  can  be  obtained  by  considering  the  integral 
curves 


(p  -  J 
o 


(i.n) 


in  the  P  -  ft  plane  (Figure  8).  As  is  to  be  expected,  the  shape  of 
the  trajectory  depends  on  whether  E  is  less  than,  equal  to,  or 
greater  than  1/32. 
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Figure  8.  Integrals  curves  in  the  P  -  f>  plane  for 
equatorial  orbits. 


For  E  >  1/52  all  trajectories  run  off  to  infinity,  and  no 
periodic  solutions  exist.  For  E  »  1/52,  the  circle  p  s  2  is 
a  periodic  orbit  (in  the  x  -  y  plane).  Moreover,  one  trajectory 
spirals  into  this  circle  from  within  and  one  from  without  (Figure  9). 
For  E  <  I/52  there  exist  two  distinct  types  of  orbits.  For 
Pq  <  P<  2,  the  orbits  are  all  periodic,  and  for  p  >2,  the  or¬ 
bits  run  off  to  infinity. 


Figure  9.  Orbits  in  the  x  -  y  plane  spiraling  into  a 
periodic  orbit,  for  £  »  1/32. 

Stormer  has  done  extensive  work  in  calculating  other  families 
of  periodic  solutions  [1*» ],  A  general  method  to  obtain  all  these 
periodic  motions  vas  presented  by  De  Vogelaere  [*♦].  We  shall  ob¬ 
tain,  as  a  corollary  to  our  result  on  the  existence  of  almost 
periodic  solutions,  infinitely  many  periodic  solutions,  vhlcb  like¬ 
wise  will  lie  in  the  oval  like  region  of  Figure  3 . 

The  behavior  of  the  equatorial  orbits  for  small  excursions  out 
of  the  equatorial  plane  may  be  examined  by  perturbation  methods. 
Expanding  V  as  a  power  scries  in  s  and  applying  Hamilton's 
equations,  one  obtains 
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S  + 


0  + 


0(2) 


(1.12) 


P+  ir(P-l)(2-P)  =  0  +  ^(z2).  (1.13) 

P? 

If  the  terms  of  second  and  higher  order  in  z  are  neglected,  the 
solutions  to  equation  (1.13)  are  the  equatorial  orbits,  and  are 
therefore  known  periodic  functions  of  time  (for  p  <  2,  E  <  1^2). 
Equation  (1.12)  then  becomes  Hill's  equation.  Its  solution  can  be 
written  in  the  form 


z  -  Ce^ft)  +  De’^V-t)  (1.14) 

where  C  and  D  are  arbitrary  constants  and  *(t)  is  periodic  in 
time  with  the  same  period  as  p(t).  The  constant  ft,  the  character, 
istic  Poincare  exponent  determines  the  stability  of  the  orbit.  It 
can  be  only  real  or  purely  imaginary.  If  n  is  real  and  4  0, 
the  motion  in  the  z-direction  grows  (within  the  approximation  made) 
without  bound.  If  ft  is  purely  imaginary  the  motion  is  bounded 
(for  time  intervals  in  which  the  neglected  terms  have  negligible 
effect)  and  is  therefore  stable.  The  behavior  of  ft  as  a  func¬ 
tion  of  has  been  studied  by  De  Vcgclaere  f5)  who  finds  that 
alJ  orbits  are  stable  for  7^  >  1.3137. 

In  Section  3,  we  shall  show  the  existence  of  a  family  of  two 
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dimensional  invariant  tori  in  the  four  dimensional  phase  space  of 
a  particle.  Any  trajectory  starting  between  two  such  tori  can 
never  escape,  and  will  remain  trapped  between  these  two  tori  for¬ 
ever.  This  is  a  much  stronger  result  than  De  Vogelaere' s.  However, 
our  result  is  only  valid  for  valu' s  of  7^  much  greater  than  1. 3137. 

2.  Almost  periodic  Motions,  Moser* s  Theorem 

(a)  We  will  be  concerned  in  this  paper  with  Hamiltonian  systems 
of  two  degrees  of  freedom,  and,  in  particular,  with  proving  the  ex¬ 
istence  of  quasi-periodic  solutions  for  systems  close  to  integrabls 
oms.*  To  this  end  we  turn  now  to  the  discussion  of  two  geo¬ 
metrical  theorems  which  will  be  basic  for  the  following.  These 
statements  refer  to  area  preserving  mappings  defined  in  an  annulus 
in  the  plane.  How  the  reduction  of  the  differential  equations  to 
a  mapping  can  be  carried  out  we  will  see  later  on. 

To  describe  an  annulus  in  the  plane  we  use  polar  coordinates 
2  2 

R  ■  x  +  y  and  the  polar  ingle  ©.  The  annulus  is  given  by 

1  s  3  f  2 


and  the  area  element  by 


There  is  a  theorem  by  Kolmogorov  and  Arnold  [2]  guaranteeing  the 
continuation  of  quasi-periodic  motions  under  perturbations 
cf  the  Hamiltonian.  However,  this  theorem  does  not  quite  apply 
to  our  case  since  the  second  frequency  will  be  small.  For 

this  singular  case  we  resort  to  Moser's  theorem. 


dxdy  =  ^dRd6. 

At  first  we  consider  a  simple  type  of  mapping  which  we  shall  call 
a  "twist  mapping" 


ex  »  0  +  r(R) 

which  clearly  preserves  the  area  element  as  well  as  concentric 
dries  R  =  constant.  Each  of  these  circles  is  rotated  by  an 
angle  ?(R),  which,  in  general,  depends  on  the  radius.  It  will  be 
a  basic  assumption  for  the  following  that  7(R)  is  not  a  constant, 
or,  more  precisely,  that 

in  1SR5  2.  (2.1) 

The  properties  of  this  mapping  are  easily  understood.  Each  circle 
R  =  constant  for  which  7/2T  ■=  p/q.  is  rational  consists  of  fixed 
points  of  M^.  Each  circle  for  which  7  and  2v  are  incommen¬ 
surable  is  densely  covered  by  the  images  of  the  iterates  M^(q  = 
1,2,,..)  of  any  point  on  this  circle. 

Our  main  objective  will  be  to  discuss  a  mapping  M  which  is 
c.  ose  to  the  twist  mapping  M  .  Therefore,  we  consider  a  mapping  M 
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R1  =  R  +  €f(R,0,c) 

(2.2) 

0-L  =  0  +  7(R)  +  eg(R,0,e) 


j 

where  f,g  are  assumed  to  have  the  period  2rr  in  0.  This  mapping 
is  defined  in  the  annulus  1  %  R  S  2  hut  need  not  map  this  annulus 
into  itself. 


Theorem  (Poincare-Birkhoff).  Let  y * (R)  4  0  and  let  he  area 
preserving  for  all  €  in  the  sense  that 


/  Rd0  =  /  Rd0+ 

C  MC 

j 

for  any  closed  curve  C.  Given  any  rational  number  p/q  between 
2rr7(l)  and  2wy(2)  there  exist  2q  fixed  points  of  satis¬ 
fying 


V* 

@  =  0  +  2irp 
< 1 


provided  €  is  sufficiently  small. 

This  theorem  is  actually  a  very  simple  version  of  the  celebrated 


+For  simplicity  we  write  M  in  place  of  M  • 
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and  much  deeper  fixed  point  theorem  by  Poincare  and  Birkhoff  in 
which  no  smallness  assumption  is  required.  The  present  theorem 
is  a  consequence  of  the  implicit  function  theorem^  One  sees 
immediately  that  the  q^*1  iterate  of  M  has  the  form 

\  =  R  +  ^(e) 

0^  =  0  +  qy(R)  +  ^(e) 

and  —  by  the  implicit  function  theorem  —  there  is  a  unique  R  = 
F(0,e)  satisfying 


0=6+  2-jrp 

q 

for  sufficiently  small  e.  This  solution  R  =  F(0,e)  represents 
a  starlike  curve  C  which  is  mapped  "radially"  by  m'1,  since  0^ 
and  0  differ  by  an  integral  multiple  of  2ir.  Since  the  mapping 

q  r 

M  and  hence  M  preserves  the  area  /  Rd0  it  follows  that  C 

C 

and  the  image  curve  M^C  intersect  in  at  least  two  points.  These 
points  are  clearly  the  desired  fixed  points.  If  P  is  one  fixed 
point,  MP, provide  q  -  1  further  fixed  points  and  the 
theorem  is  proven. 

It  is  important  to  observe  that  we  had  a  continuum  of  fixed 
points  of  for  e  «  0  while  for  e  >  0  we  predict  only  a 
finite  number  of  fixed  points.  In  fact,  it  can  be  shown  by  examples 
that  in  general  the  curve  of  fixed  points  breaks  up  into  a  finite  set 


of  fixed  points 


25 


Now  we  turn  to  the  question  of  what  happens  under  perturbation 
to  those  circles  R  =  constant  for  which  y(R)/2ir  is  irrational. 
This  is  the  content  of  the  following  theorem  which  will  not  he 
proven  here. 

Theorem  (Moser  [ll]).  Let  f  (R)  ^  0  and  let  any  curve  C  sur¬ 
rounding  R  =  1  and  its  image  curve  MC  intersect  each  other. 

The  function  f,g  are  assumed  to  be  sufficiently  often  differ¬ 
entiable.  Then,  for  sufficiently  small  e  there  exists  an 
invariant  curve  T  surrounding  R  =  1.  More  precisely,  given 
any  number  o  between  y(l)  and  y( 2)  incommensurable  with  2tt, 
and  satisfying  the  inequalities 


|£- 
1  2tt 


-  fi  s  ourJ/2 


for  all  integers  p,q,  and  some  constant  c  >  0,  there  exists  a 
differentiable  closed  curve 


R  *  F(*,e) 

6  a  o  +  ,e) 


(2.5) 


with  F,G  of  period  2  tr  in  0,  which  is  invariant  under  the 
mapping  Mg  --  provided  e  is  sufficiently  small.  The  image 
point  of  a  point  on  the  curve  (2.3)  is  obtained  by  replacing 
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0  by  <t>  +  a). 

In  actual  fact  this  theorem  provides  not  only  one  but  infinitely 

many  invariant  curves.  We  observe  that  the  intersection  property 

for  curves  mentioned  in  the  above  theorems  is  certainly  satisfied  if 

the  area  /  Rd0  is  preserved.  Namely,  if  C  is  any  curve  surround- 
C 

ing  R  =  1  so  will  the  image  curve  surround  R  =  1.  If  the  curves 
C  and  MC  would  not  intersect  each  other  one  of  them  would  lie  in¬ 
side  the  other  and  the  areas  /  Rd@  and  /  Rd0  could  not  agree. 

C  MC 

For  later  applications  we  need  an  extension  of  the  previous 
theorem!  If  the  mapping  (2.2)  is  replaced  by 


Rx  =  R  +  e°f(R,0,e) 

e1  *  0  +  a  +  epy(R)  +  €°g(R,0,e) 


(2.4) 


where  OS  pea,  then  the  conclusions  of  the  previous  theorems  re¬ 
main  true.  The  essential  point  is  that  the  perturbation  term  is 
small  compared  to  the  "twist"  eP7(R). 

To  apply  Moser’ s  theorem  to  prove  the  existence  of  invariant 
tori  for  a  Hamiltonian  system  of  two  degrees  of  freedom,  we  first 
approximate  the  Hamiltonian  H,  (if  possible)  by  an  integrable 
Hamiltonian  HQJ  i.e.  we  write 

H  -  Hc(R1,R2,£)  +  €iH1(R1,01,R2,92,€) 


(2.5) 


where  o>,  =  ciK^cR^  is  of  order  one,  and  the  frequency  ratio 
ci^/u^  varies  over  a  region  of  order  e  with  k  <  i(o^  =  dHo/&g). 
On  the  energy  surface  H  =  c  we  solve  for 


(2.6) 


Using 

*2  =  R, 


0^  as  independent  variable  instead  of  t  and  setting 
=  we  find  *ram  Hamilton' s  equations  that 


dR  "H0  .  d0  ^ 

®I"V 


(2.7) 


One  verifies  easily  that  on  H  =  c  these  equations  take  the  form 


dR  .  .  d0  . 

I  ■  CS  <P  .  ■  s 

d01  0  *  d©1  R 


(2.7)* 


where  <!>  is  defined  in  (2.6). 

Hie  system  (2.7)  is  again  Hamiltonian,  of  one  degree  of  freedom, 
but  non- autonomous.  To  eliminate  the  independent  variable  we  fol¬ 
low  the  solutions  from  0.^  ■  0  to  the  next  intersection  with  0^  ■ 
2tt,  This  defines  a  mapping  which  —  by  Liouville's  Theorem  — 
preserves  the  area  /  Rd0.  This  mapping  will  have  the  form 


R(&r)  -  R(0)  +  0(t*) 


0(2»r)  -  0(0)  +  a  +  cky(R)  +  ^(«*) 


(2.8) 
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with  y'  (R)  4  0.  Note  that  the  condition  y'  (R)  4  0  means  that 

the  frequency  ratio  a>  ^  varies  on  the  energy  surface  H  =  c. 

The  important  point  is  that  the  frequency  ratio  varies  over  a 

4  k 

region  which  is  large  compared  to  e  ,  i.e.  the  twist  €  7(R)  is 
small  compared  to  the  neglected  terms.  Moser's  Theorem  guarantees 
infinitely  many  invariant  curves  of  the  mapping  (2.8),  for 
sufficiently  small  e.  Any  curve  invariant  under  this  mapping 
generates  an  invariant  torus  if  we  take  all  solutions  which  issue 
forth  from  the  invariant  curve,  and  on  these  tori,  the  motion  is 
quasi- periodic  with  two  frequencies. 

In  this  case  a  quasi-periodic  motion  will  densely  cover  a  two 
dimensional  invariant  torus  in  the  four  dimensional  phase  space  of 
the  particle.  Any  trajectory  starting  between  two  such  tori  must 
always  remain  between  these  two  tori.  (Warning:  This  is  not  true 
for  n  >  2.)  This  provides  a  powerful  tool  for  proving  the  sta¬ 
bility  of  periodic  orbits.  In  the  following  section  we  shall  show 
the  existence  of  quasi-periodic  motions  in  the  earth's  magnetic 
field  which  lie  near  the  equator.  All  trajectories  near  the 
equator  which  are  caught  between  two  such  tori  must  remain  near 
the  equator  forever.  This  is  the  stability  result  for  "near 
equatorial"  orbits  which  we  mentioned  in  the  preceding  section. 

(b)  We  would  now  like  to  show  how  one  goes  about  approximating 
H  to  higher  and  higher  order  by  an  integrable  Hamiltonian.  The 
method  we  shall  present  is  known  as  the  lindstedt  method  [13). 
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For  simplicity  we  shall  restrict  ourselves  to  systems  with  two 
degrees  of  freedom. 

Consider  a  Hamiltonian  H  of  the  form 

H1  =  Vl  +  V2  +  €Vxi>X2»yl'y2)  +  +  (2.9) 

where  x^  is  conjugate  to  yi  and  H  has  period  2ir  in  y 
and  y2,  i.e. 

H(x,yi+2rr)  =  H(x,yjL). 

We  would  like  to  find  new  canonical  variables  x£#yj,  so  that  H 
will  be  independent  of  the  angular  variables  y^,y£  through  terns 
of  order  n.  To  this  end  we  consider  a  generating  function  of  the 
fora 

s  ■  yixi  +  y2*2  +  €Si(xi»xa»yi»y2^  +  *  •••  (2*10) 

with 

xi  "  »  yi  ■  •  (2.io)' 

If  we  denote  the  Hamiltonian  H  expressed  in  term  of  the  primed 

variables  by 
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H(x’,y’) 


+  Ui^'2  +  cR^xy  +  ^2  + 


(2.11) 


then  from  (2.10) 


x' 

X2 


*1 


yl>y2) 


-fl. 


(2.12) 


Equating  terms  of  order  €  in  (2.12)  we  see  that 


cB  at 

(^1  +  o>2  “  ^i(xi»x2^  ~  ^1  (X1>X2* *  (2*15) 


In  order  for  S.^  to  he  periodic  in  y1  and  y2  we  must  require 
that  the  right  hand  side  of  (2.13)  have  mean  value  zero.  Hence 

2 ir  &r 

Wa*  '  7^5  II  Hi(xi'^yi>y2)dyidy2>  (**1M 

i.e.  JtL  is  the  mean  value  of  H1»  We  still  cannot  solve  for  Sx 
since  a  Fourier  series  expansion  will  contain  the  small  divisors 
^1%  +  therefore  require  that  the  frequencies  satisfy 

the  Infinitely  many  Inequalities 


l(J,®)l  *  rlJl*’ 


(s.») 


for  all  Integers  with  |  j|  «  jjJ  ♦  |jg|  >  0,  and  with  same 

constant?  y  and  t  >  l.  If  the  right  hand  side  of  (2,13)  has 
the  Fourier  expansion 
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.  r. 

kio 


V*') 


ei(*,y) 


then  is  given  by 


S 


1 


V 

i— 

k^o 


\(x,) 

TTJvSI 


ei(^,y)  . 


To  determine  we  note  that  equating  terms  of  order  k  in 
(2.12)  yields 


a&  as 


(2.16) 


where  depends  only  on  0  >  i  *»  k  -  1,  Thus  is 

determined  by  setting  the  mean  value  of  the  right  hand  side  of 

(2.1b)  equal  to  xero,  and  then  is  determined.  To  eliminate 

the  y  dependence  of  H  through  order  n,  we  simply  truncate  the 

series  for  S  after  «nS  . 

n 

The  above  method  may  bo  extended  to  the  ease  where  the  lowest 
order  term  H  of  H  is  given  by 


Hc  -  VX1'X2*' 


c 


We  simply  choose  x*  ■  x  so  that 
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are  rationally  independent  numbers  satisfying  (2.15),  and  then  re¬ 
strict  the  variables  x’  to  a  small  neighborhood 

Finally,  if  the  frequency  o>2  is  zero  to  lowest  order,  we 
first  eliminate  the  y^  dependence  of  H  through  order  n.  Then, 
if  we  can  find  new  variables  so  that  is  also  independent  of 
y^,  we  may  eliminate  the  y?  dependence  through  order  n. 

3.  Main  Result.  Existence  of  Quasi- Periodic  Motions 

(a)  In  this  section  we  will  prove  the  existence  of  quasi- 
per-'odic  solutions  of  the  Storaer  problem.  These  motions  will  all 
lie  in  the  oval  like  region  of  Figure  5  and  will  satisfy  H  «  1/32. 
Moreover,  these  motions  may  penetrate  arbitrarily  close  to  the 
dipole. 

(b)  Dipolar  coordinates.  Since  our  intuitive  idea  of  the 
motion  is  a  gyration  about  a  line  of  force  and  an  oscillation  a- 
long  the  line  of  force,  it  is  natural  to  introduce  new  " orthogonal" 
coordinates  a(o,z),  b(p,z)  such  that  the  curve  a(p,z)  «  const, 
defines  a  magnetic  line  of  force.  The  magnetic  field  lines  are 
given  by  Uie  equation 

a  -  ,  (a  •  const.). 

cos^X 


(3.1) 
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From  the  equation 


da  3b  da  3b  A 

stTp  +  ss*  0 


(3.2) 


we  find  that 


»(p,*>  .  ^  .  (3.3) 

r 

(Any  function  of  b  together  with  a  provide  an  orthogonal  co¬ 
ordinate  system. )  The  new  canonical  variables  may  be  obtained 
from  the  generating  function 


F(P>z*Pa>Pb)  »  *(p,z)pa+  (3.U) 

by  employing  the  standard  relations 


<F  da  db 

•- ■3T  ' 

& 


(3.5) 


In  terms  of  these  new  variables  the  Hamiltonian  H  of  (1.9)  takes 


the  form 


where 


2  _  cos^A  ,2  _  a^cos^^A 
a  LOsin^A  ’  °  l+3sin*”A  * 


(3.7) 


It  is  understood  that  sin  A  and  cos  A  are  to  be  expressed  in 
terms  of  a  and  b. 

We  now  wish  to  restrict  ourselves  to  a  region  in  phase  space 
where  the  energy  H  will  be  small.  This  is  to  conform  with  our 
notion  that  in  the  course  of  one  gyration  about  its  guiding  field 
line,  the  particle  should  see  an  approximately  constant  magnetic 
field.  Thus  we  consider  the  change  of  variables 

2 

a  -  1  =  e  a  b=ep 

2  5  »  (5’8) 

\  -  5  pa  -  €  pe 

2 

where  e  is  a  small  parameter.  The  condition  a  -  1  =*  e  a  means 
that  we  require  the  particle  to  remain  near  the  guiding  field  line 
r  =  cos  2A.  Our  system  will  remain  Hamiltonian  if  we  take 

=  H/'€^  (3.9) 

Our  next  step  is  to  try  and  approximate  H  by  an  integrable 
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Hamiltonian.  To  this  end  we  first  show  that  X(a,b)  is  an 
analytic  function  of  a  and  b  for  b  small.  Squaring  equation 
(3.1)  and  multiplying  by  (3.3)  yields 


sin  X 

(l-sin2*)2 


(3.10) 


The  derivative  of  the  right  hand  side  of  equation  (3.10)  with  re¬ 
spect  to  X  is  one  at  X  =  0.  Therefore,  we  are  guaranteed  that 
X(a,b)  is  an  analytic  function  of  z  =  a^b  for  |  z|  small. 
Hence,  we  may  write 


p  2  2  4  2  4 

sin  A  =  €  p  +  !♦€  (a£  -0  )  +  G^a,^) 

(3.11) 

cos'^X  =  1  +  3e2p2  +  6e1+(aa02-p2+)  +  G2(a,p,€) 


where  G.  can  be  written  in  the  f  €  G  •  (a,P,e)  with  C. 

i  i  '  x 

analytic  in  the  variables  a,P,e.  These  expansions  are  trivially 
derived  from  the  relation 


(3.12) 


The  Hamiltonian  (3.9)  nay  now  be  written  in  the  form 


•■VVVH( 


(3.15) 


where 


3*4 


2  2 

a  +p„ 


H2  »  ^(6pSp^-l.a3+Pp+3eV) 

Ht  =  ^-(^%2-XkPQ-2aPp+?e%p+2a-2a%k) 


and. 

Hg  =  €6H6(a,P,Po,Pp,0 

with  analytic  in  all  its  variables. 

We  shall  now  show  that  Hq  +  can  be  transformed  into 

g 

an  integrable  Hamiltonian,  modulo  terms  of  order  €  .  Firstly, 
define  new  coordinates  R, 0  via  the  formula 

a  ■  V^R  sin0,  pq  =  V2R  cos0.  (3.14) 

From  their  definition,  R  and  0  are  canonical  coordinates.  The 

k 

magnetic  moment  M  of  the  particle  will  be  proportional  to  e  R 
(to  be  shown  in  Section  4),  Note  also  that  HQ  =  R>  and  R  is 

p 

constant  to  order  e  .  Next,  we  employ  Lindstedt' s  method  to 
average  out  the  0  dependence  of  H  to  order  e6,  i.e.  we  define 
new  variables  R'>8',0',Pg  so  that  H  is  independent  of  0 
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through  order  The  frequency  u>2  of  (2.9)  is  zero,  while 

<0,  =  1.  In  the  notation  of  Section  2, 


l  2fT  2  ^  p 

H2  =  ^Jr  /  H2d0  "  V9R,(P’>  +  1 

o 


(3.15) 


and 


i  2rr  5  2  ^Hp  2  35p  ^Bp 

\  =  ^!  [\+e  *5T-sr+e  ^^r140*  (3*l6) 

o  p  p 


The  generating  function  S0  =  32(R',e,P,p£)  is  given  by 


=  -  sin  20  -  4(2R* )3/2(-cos  9  +  ^^)].  (3.17) 


Performing  the  integration  in  (3.16)  we  find  that  the  Hamiltonian 
(3.13)  may  be  written  in  the  form 


Ho  +  H2  +  HU  +  H6 


(3.18) 


where 


H  -  R 
o 
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(we  have  suppresses  the  primes  for  convenience),  and  Hg  =  ), 

For  fixed  R  the  curves 


2  2 

QRfJ  +  pQ  =  const. 
P 


are  ellipses  in  the  g,  pc  plane.  These  curves  may  he  transformed 

p 

into  circles  (with  the  same  area)  by  the  generating  function 


F(e,p,Rpp£)  =  (9Rl) 


1/4, 


Bp  + 


0R. 


(3.19) 


Setting 


P*  =  \/5^  sin  Sg,  p»  =  cos  @2  (3.80) 


we  see  that  is  independent  of  9^.  Hence,  we  may  employ  the 

Linds tedt  method  to  average  out  the  9^  dependence  of  H  to 
order  e^.  In  terms  of  new  canonical  variables  which  we  again 
call  Hamiltonian  H  now  assumes  the  form 


o  n  €k  2  13R2 

H  =  R^  +  3s  RgvR^  +  ’TT^”’^’R1  +  *"'g"  *)  + 


(3.21) 


where  Hg  ■  e^?g(R1,Rg,01,02,€),  with  !g  analytic  in  all  its 
variables.  Thus,  we  have  succeeded  in  approximating  H  by  an 
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integrable  Hamiltonian 


F  -  «  -  H6 

to  order  e* .  We  now  apply  Moser's  Theorem  to  prove  the  continua¬ 
tion  of  quasi-per iodic  motions  under  the  perturbation  Hg.  As 
described  in  Section  2,  we  solve  for  R^  =  e)  on  the 

energy  surface  H  =  c,  and  take  0^  instead  of  t  as  independent 
variable.  We  then  follow  the  solutions  from  0^  =  0  to  their 
next  intersection  with  0^  =  2rr.  This  defines  an  area  preserving 
mapping  which  we  denote  by  M.  Since 

33“  "  ’  |r  **  R2  +  (5*22)  • 


and 


bV2.c1/2 


2  + 


(3.22)' 


the  mapping  M  has  the  form,  in  the  coordinates  R  *  R^,  0  »  0^ 


R  m  R(2rr)  =  R  +  e  f(R,0,e) 


m: 


0  *  0(2ir)  -  0  +  3e2c1/2  -  ^  e^R  +  «6g(R ,0,e) 


(3.23) 


where  R  =  R(0),  0  -  0( 0)  and  f,g  are  analytic  in  the  variables 
R,0,€.  This  mapping  is  of  the  form  (2.4)  with  p  =  4,  o  =  6,  and 


r  00  =  -  22  +  o. 


Hence,  Moser* s  Theorem  applies. 

Thus,  we  have  established  the  existence  of  quasi-per iodic 
motions  with  two  frequencies  for  sufficiently  small  e.  Since  we 
restricted  b  to  be  small,  all  these  motions  must  lie  near  the 
equatorial  plane.  Moreover,  these  motions  gyrate  tightly  about 
the  guiding  field  line  r  »  cos2^  (or  r  *  r^cos2*  in  our  old 
coordinates).  A  typical  motion  is  illustrated  in  Figure  10. 


Figure  10.  A  quasi-periodic  motion  in  the  p  -  z  plane 
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These  orbits  all  cross  the  magnetic  field  line  at  an  an^le  near 

90°  since  the  velocity  of  the  particle  parallel  to  the  magnetic 

field  is  much  smaller  than  the  total  velocity.  This  follows  from 

the  fact  that  p^  =  e^p^  is  of  third  order  small  while  p&  = 

2 

e  d  is  only  of  second  order. 

*a 

(c)  Quasi-per iodic  motions  which  penetrate  arbitrarily  close 
to  the  dipole. 

Our  goal  now  is  to  prove  the  existence  of  quasi-periodic 
motions  which  need  not  lie  near  the  equatorial  plane.  To  this  end 
we  consider,  instead  of  (3.8)  the  change  of  coordinates 


a  -  I*  ea  J  b  =  tP 

p»  ■  %  ;  v  epP 


(3.2*0 


where  again  €  denotes  a  small  parameter.  Let  M  and  N  be  two 
fixed  constants,  with  N  very  large.  We  then  consider  all  vari¬ 
ables  to  be  complex,  and  restrict  ourselves  to  the  region  T 
defined  by 


M  +  |Pal  +  |ppl  *  M 

| Re  b|  *  N 


(3.25) 


|Im  b|  s  S(M,N,e) 


bO 

where  5  depends  on  M,N,  and  e  and  will  be  chosen  so  that  the 
Hamiltonian  (3.6)  is  analytic  in  the  variables  a,b,pa,pt  in  the 
region  T,  for  €  sufficiently  small.  To  find  the  singularities 
of  H  in  the  complex  4-space  we  consider  again  the  equation  (3.10) 


sin  X 
(l-sin2*)2 


With 


and  y 


sin  A(a,b),  we  observe  that 


dz 

3y 


i+3y2 

(1-yV 


(3.26) 


Hence  the  only  singular  points  are  y  =  tl,  y  =  t  •  It  is 

clear  that  for  fixed  M  and  N,  6  may  be  chosen  so  that  sin  X  J 
tl  for  |  La  bj  5  6,  and  «  sufficiently  small,  ltie  points  y  ■ 

t  i'y/^  correspond  to  the  points 


z  ■  a^b  «  ti 


Letting 


2 

a  ■ 


ai  + 


b  »  b^  +  ibg, 


we  see  that 


a^b  ■  a^  -  Sgbg  +  *(aib2+a2bl^ 


(3.2?) 


For  fixed 


and  €  still  further 


4l 

*•  and  N,  we  now  restrict  6 

so  that 


lalVa2bll  “I! 

thus  excluding  the  points  z  =  t  new  Hamiltonian 


p  p 

H(a,Pa,P,Pp,€)  =  H(a,h,Pa,Pb)  =  ^(-g  +  ^|)  + 


P2  P2  2 

:  l  Vt  (3*aB) 

h^  h^  2(l+€a)4cosbX 


is  now  analytic  in  sill  variables  in  the  region  T.  Note  that  al¬ 
though  N  is  fixed,  it  may  be  chosen  as  large  as  desired.  The 
points  (a  «  1  +  eor,  b  ■  all  lie  very  close  to  the  dipole,  and 
approach  the  dipole  as  N  -*  *. 

(d)  Our  next  task  is  to  Approximate  the  Hamiltonian  by  an  in- 
tegrable  one.  We  cannot  expand  H  in  powers  of  so  and  d?  as 
we  did  previously,  since  now  :•  "•  l'e.  Instead  we  expand  H  in 
rover 8  of  (a-1)  «  so  only.  Since  sin  X  and  cos”1  are  analytic 
in  T  we  may  write 

1+3  sin2X  .  KLCb)  +  F^b) 
cos”bX  ■  C^(b)  +  Fg^b) 


(3.29) 


where  Fi(a,b)  =  <^((a-l)).  (^((a-l)X)  denotes  an  analytic 

function  which  vanishes,  together  with  its  first  k  -  1  derivatives 
with  respect  to  a,  at  a  =  1. )  The  H?*miltonian  (3.28)  may  thus 
be  written  in  the  form 


H  =  H  +  H. 
o  ± 


(3.30) 


where 


Ho  ■  7<KiVa2> 


■e 


(3.30)' 


and  with  Ux  analytic  in  all  its  variables. 

The  arguments  of  the  functions  Cx  and  are  b  «  c0. 

(e)  The  next  logical  step  is  to  transform  the  Hamiltonian  Hq 
into  an  integrable  one,  modulo  terms  of  order  c.  We  start  by 

2 

finding  new  variables  a',p^  80  that  H0  8  faction  of  (ft* )  + 
(P»)2  alone.  These  variables  may  be  obtained  from  the  generating 

v’0f 

function 

F(o.e.^Pa)  •  k'1  ^  (3-W) 

where  the  argvents  of  the  functions  Kx  and  f  are  eft.  It 
would  suffice  (for  the  purpose  stated  above)  to  take  f  ■  «ft.  How. 
ever,  a  Judicious  choice  of  f  will  enable  us  to  express  HQ 


^3 


explicitly,  and  thereby  simplify  most  of  the  later  calculations. 
The  new  canonical  variables  a' ,  p^ £ ' , are  determined  from  the 
relations 


a?  K-i'4 ,  .  flt  a?  _  f 

75  “  K1  Pa  >  P  =  ^  -  « 


a  =  W =  K1  a 


fr1  .(1)  ,  ev-l/Uv(l)  , 

-  -JT  s  f  -  lft.  K’  OT 


(3.32) 


1  % 


where  and  f^  denote  the  derivatives  of  the  functions 

and  f  with  respect  to  v  »  €0. 

A  point  in  the  p  -  z  pxane  is  determined  from  the  coordinates 
a’ in  the  following  manner;  Given  a*  and  p*,  the  point  lies 
on  the  magnetic  field  line  r  *=  a  cos^>,  where 

a  *  1  +  o*n+3(€e')2]I/U.  (3.33) 

the  coordinate  b  on  this  magnetic  field  line  is  then  found  from 
the  relation 


f(b)  -  f(«f)  -  €?•.  (3.3M 

We  cannot  solve  (3. 3t)  exactly  for  b.  However,  for  the  f  we 
shall  determine  it  will  be  seen  that  as  cf  ’  increases  (decreases) 
to  ♦!(-!),  b  increases  (decreases)  to  ♦•(-•)*  and  that  the 
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latitude  \  for  a  =  1,  is  precisely 

sin"1(e^' ). 

Thus,  p* (t)  essentially  describes  the  motion  along  a  line  of 
force,  for  the  coordinate  a  near  one.  In  terms  of  the  primed 
variables  the  Hamiltonian  (3.30)  assumes  the  form 


H  = 


(P^)2+  C^ff')^)2]  +  Hx  , 


(3.35) 


with  H.  =  <^(e). 

(f)  We  shall  now  show  that  the  Hamiltonian  H  is  an  in- 
tegrable  Hamiltonian,  i.e.  we  shall  find  new  canonical  coordinates 
R,6,J,(t>  so  that  Hq  =  Hq(R,J).  To  this  end  we  define  canonical 
coordinates  R,0  in  place  of  a'  ,p^  by  the  relations 

a'  =  y/2R  sine,  p'  =  V®  cose,  (3.36) 

CL 

Hq  is  now  independent  of  9.  Hence  the  coordinate  R,  which  we 
shall  later  show  (see  Section  4)  is  proportional  to  M/e  ,  is  a 
constant  of  the  motion  to  lowest  order  in  e.  To  lowest  order,  0 
is  given  by  the  equation 


t 


(3.37) 
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But  the  quantity  on  the  right  hand  side  of  (3.37)  is  the  magnitude 
of  the  magnetic  field  along  the  line  of  force  r  =  cosSv  (see  (1.1), 
(1.2)  and  (3.29)).  Thus,  we  have  confirmed  our  intuitive  notion 
that  to  lowest  order  the  particle  gyrates  about  its  guiding  field 
line  with  the  cyclotron  frequency  oj^. 

It  is  convenient  at  this  point  to  introduce  a  new  time  scale  t 
so  that  the  frequency  becomes  one,  i.e.  the  particle  gyrates 
about  its  guiding  field  line  in  the  new  time  scale  with  constant 
frequency.  This  is  easily  accomplished  by  setting 


t  =  f  co,dt. 


(3.38) 


Our  trajectories  are  now  the  zero  energy  solutions  of  the  Hamiltonian 


F  =  i-(H-h%  =  F  +  F 

/  in  4 


U) 


(3.39) 


where  h  is  ^he  constant  value  of  the  Hamiltonian  (3.35)  and 


1  /  2,  »  2,  ,\2 


F0  =  R  +  (pp) 


2h) 


H1 

F.  =  —  . 
1 


(3.39)' 


Thus,  to  lowest  order  we  have  decoupled  the  motion  in  or’,p^  and  in 

p*  p*.  This  puts  into  evidence  the  integrability  of  the  truncated 
P 

system. 
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In  this  approximation  the  motion  in  the  P' ,P^  plane  is 


described,  by  the  Hamiltonian 


*  =  K(f,)2(pp)2-l[ 


(5.4o) 


which  we  now  calculate.  For  this  purpose  we  observe  that  (a±  = 
c  ^  and  f*  all  are  expressed  in  terms  of  f(b)  =  f(cP)  =  eP’ 
if  we  define  f (b)  by 


(1  -fT 


(3.4l) 


Therefore,  <D  can  be  expressed  in  explicit  algebraic  form.  To 


show  this  write 


sin  A  «=  f^(b)  +  ^((a-1)) 


where  we  have  f  «  f  (see  a^b  =  '  'g  — g) »  Then  compute 


(l..sin  a) 


a,  =  ¥-■  with  7  =  eP’ 

1  (1-r  r 


„  .  «  U-r2)? 


Ojj  Vl+372 


I 
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Hence 


svG3?  1+3?‘ 


(3.42) 


The  motion  in  the  £'  pi  plane  derived  from  the  Hamiltonian 

p 

Fq  is  qualitatively  determined  by  considering  the  level  curves 
(see  Figure  11) 


-  -  -‘'(I-?2)3  .  1  (1-72)3 


(pA)‘ 


in  the  pi  plane.  For  c  >  0  these  curves  are  not  closed  and 

'  p 

asymptotically  approach  the  lines  y  =  ±1.  The  curves  corresponding 
to  c  =  0  are  hyperbolas  defined  by 


IP  P 

T;(Pg)  -  3h7  --  h 


(3.44) 


U8 


Figure  11.  Level  curves  in  the  y,  p£  plane  of  <J> 


=  const. 


For  c  <  0,  the  components  of  these  curves  in  | y\  <  1  are  all 
closed,  and  lie  inside  the  domain  bounded  by  y  =  tl  and  the  hyper¬ 
bole,  with  the  point  s  in  Figure  11  determined  from  the  equation 


h  _  (l+3s 

c  ~ 


2^1/2 

r 


(1-s  ) 


OM) 


The  quasi-periodic  orbits  obtained  previously  correspond  to  small 

oscillations  in  the  P',p'  plane  about  the  equilibrium  correspond- 

p 

ing  to  equatorial  orbits.  Since  we  are  only  considering  the  zero 
energy  solutions  of  F,  the  constant  c  is  essentially  -R.  Hence 
P' (t)  and  Pp(T)  will  be  periodic  functions  of  t  if  we  neglect 
the  term  F1  in  (3.39).  Thus  we  have  .justified  our  intuitive  idea 
that  the  guiding  center  of  a  particle  oscillates  along  a  line  of 
force  between  two  mirror  points.  This  motion  is  indicated  sche¬ 
matically  in  Figure  12.  To  lowest  order,  the  particle  mirrors  at 
a  latitude  A  =  sin-1s  where 


Figure  12.  Schematic  description  of  guiding  center  motion. 

Thus,  the  smaller  R,  the  closer  the  particle  approaches  the 
dipole.  It  is  also  clear  from  Figure  11  and  the  discussion  con¬ 
cerning  the  coordinates  a' that  there  exist  orbits  of  the 
unperturbed  system  Fq  in  the  region  |b|  s  N  which  achieve 
|  b|  >  3N/4. 

Finally,  since  the  curves  (3.43)  are  closed  for  c  <  0,  we  may 
introduce  the  familiar  action-angle  variables  J,  <t>  in  place  of 
P'  p',  where  the  action  J  is  given  by 

p 


ky/2 


fh(l+3x  )  +  0 


^  .  „q*5xg)3/g|1/g 


(1-1O7 


cbc. 


(3.1*7) 
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(J  is  simply  the  area  of  the  closed  curves  (3.43)  in  the  0',p^ 
plane.)  In  terms  of  the  variables  11  the  Hamiltonian  Fq 

takes  the  form 

Fo=R  +  c(eJ).  (3.48) 

Thus,  we  have  succeeded  in  showing  that  the  Hamiltonian  F  is 
’  o 

integrable.  The  two  frequencies  of  motion  determined  by  Fq  are 

o 

“i  =  ir 

dF 

o 

“s  ~  IT 

where  c'  denotes  the  derivative  of  c  with  respect  to  eJ.  In¬ 
cidentally,  the  fact  that  m0  is  proportional  to  e  is  in 
agreement  with  our  notion  that  the  particle  oscillates  slowly  a- 
long  a  line  of  force,  while  gyrating  rapidly  about  it.  Note  also 

that  the  action  J,  which  is  the  longitudinal  adiabatic  invariant 
1 a  2 

divided  by  e  ,  is  constant  t.o  lowest  order. 

We  now  apply  Moser' s  Theorem  to  prove  the  existence  of  quasi- 
periodic  motions  for  the  full  Hamiltonian 

F  =  R  +  c(eJ)  +  eF^R,  eJ,a,V)  (5.50) 

where  F^  is  analytic  in  all  its  variable.  On  the  energy  surface 


=  1 

(3.49) 

=  ec' 
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F  =  0  we  solve  for  R  =  0(cJ, 0,  $,e)  and  take  0  insteau  of  t 
as  independent  variable.  We  then  follow  the  solutions  from  0=0 
to  their  next  intersection  with  0  =  2ir.  This  defines  a  mapping 
M,  which  in  the  coordinates  r  =  eJ, $  is  given  by 

rx  =  r(2ir)  =  r  +  €2f(r,4>,e) 

m:  p  (5.51) 

<t>1  -  (t>(2r)  =  $  +  ec'  (r)  +  e  g(r, <t>,e) 

with  f, g  analytic  in  all  its  arguments.  The  intersection  property 
of  a  curve  with  its  image  curve  still  holds,  trivially.  Hence,  the 
above  mapping  M  is  of  the  form  (2.4)  with  p  =  1,  a  =  2  and 

7(r)  =  c* (r).  (3.52) 

To  apply  Moser' s  Theorem  now,  we  need  only  verify  the  condition 

7' (r)  =  c" (r)  ^  0  (3.55) 

i.e,  7(r)  should  be  a  monotonic  function  of  r.  Unfortunately, 
we  cannot  calculate  7(r)  explicitly  and  must  be  content  with 
numerical  calculations.  Figure  13  below  shows  the  graph  of  7(r) 
versus  c.  It  is  clear  that  7(r)  has  a  single  stationary  point 
at  c  =  -.70.  At  this  point  the  non- degeneracy  condition  (5.55) 


is  violated, 


Figure  13,  Graph  of  7(r)  versus  c. 

Thus,  with  the  exception  of  one  point,  Moser's  Theorem  guarantees 
the  existence  of  quasi-periodic  motions,  close  to  the  unperturbed 
motion  defined  by  the  Hamiltonian  F  ,  for  sufficiently  small  e. 

Since  some  of  the  unperturbed  motions  achieved  |b|  >  3N'U,  then  for 
sufficiently  small  e  there  will  exist  quasi-periodic  motions  for  the 
full  Hamiltonian  which  achieve  jbj  >  N/2.  By  choosing  N  aritrarlly 
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large  these  quasi-per iodic  motions  can  penetrate  arbitrarily  close 
to  the  dipole.  Of  course,  the  larger  we  choose  N,  the  smaller  we 
must  take  e,  i.e.  the  tighter  the  particle  gyrates  about  its 
guiding  field  line. 

We  now  wish  to  schematically  describe  these  quasi-per iodic 
motions  in  coordinate  space,  rather  than  phase  space.  In  the  p  - 
z  plane,  the  particle  is  confined  to  the  oval  like  region  in 
Figure  5,  with  V  very  small,  and  it  rotates  slowly  around  the 
z-axis  in  accordance  with  the  equation 

•  (t)  -  /  (i  -  £r)dt’  +  const.  (3.5M 

p  td 

It  will  be  shown  in  the  following  section  that  the  quantity 
p^h^  in  (3.6)  is  v„  ,  where  v„  is  the  velocity  of  the  part¬ 
icle  parallel  to  the  magnetic  field.  To  lowest  order  in  c  we 
then  have 


X .  -i .  $ 

v*  R  fH-3(eP»)2'! 

Fcr  R  small,  the  particle  crosses  the  equatorial  plane  (P*  »  0) 

at  an  angle  near  90°.  As  ep'  increases  to  its  maximum  value  s, 
2  2 

the  quantity  v„ /Vj_  decreases  monctonically  to  zero.  When  the 

2.  • 

particle  crosses  the  field  line  r  ■  cos  i  ■  0,  and  hence  the 
velocity  vector  lies  in  the  meridian  (0  -  z)  plane.  Thus,  the 
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particle  crosses  the  guiding  field  line  in  a  sequence  of  angles 

o 

which  monotonically  approach  90  .  This  situation  is  illustrated 
in  Figure  l4. 


Figure  l1*. 


As  the  particle  moves  tack  toward  the  equator,  the  crossing  angles 
tcgin  to  decrease  monotonically,  until  the  particle  is  moving  nearly 
parallel  to  the  field  at  the  equator.  Note  that  as  the  particle 
moves  into  regions  of  high  latitude,  it  begins  to  gyrate  rapidly, 
since  the  cyclotron  frequency  becomes  large. 

(g)  As  a  corollary  to  the  above  result  on  the  existence  of  quaei- 
periodic  actions,  we  are  automatically  guaranteed  infinitely  many 
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periodic  solutions.  We  outline  the  proof  briefly.  Each  quasi- 
periodic  motion  densely  covers  a  two  dimensional  invariant  torus 
in  the  four  dimensional  phase  space  of  the  particle.  Moreover, 
this  torus  must  lie  on  the  energy  surface  H  =  const.  We  now  cut 

i 

the  torus  with  a  surface  of  section  S,  and  consider  the  mapping  i 

P  -»M(F)  induced  on  S  by  the  differential  equations^  i.e.  M(P) 
is  the  point  where  the  trajectory  beginning  at  P  returns  to  S. 

From  the  form  of  the  Hamiltonian  F  of  (3.53),  it  is  possible  to 
introduce  coordinates  R,0  on  S  so  that  the  mapping  M  has  the 
form 


R^  =  R  +  e2f(R,6,  e) 

m: 

2 

S-L  =  8  +  ey(R)  +  e  g(r,0,e). 

A  mapping  M  of  this  form  is  referred  to  as  a  "twist"  mapping,  and 
the  Poincare-Birkhoff  fixed  point  theorem  (Birkhoff  [3])  guarantees 
the  existence  of  at  least  two  fixed  points.  Since  each  fixed  point 
represents  a  periodic  solution,  we  have  established  the  existence  of 
infinitely  many  periodic  solutions  to  the  Stormer  problem  (which  al¬ 
so  lie  in  the  oval-like  region  5,  with  H  very  small). 

4.  Quasi- Periodic  Motions  in  a  Rotationally  Symmetric  Magnetic 
"Mirror"  Field. 

(a)  We  consider  now  a  rotationally  symmetric  magnetic  field 
B,  i.e.  a  field  which  is  independent  of  the  azimuthal  angle  <t>. 


» 
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The  motion  of  a  charged  particle  in  such  a  field  can  he  described 
by  a  Hamiltonian  of  the  form 

H(p,z,pp,pz)  =  |(pp  +  p^)  +  ^  -  A(p,z))2,  (4.1) 

where  we  have  normalized  the  mass  and  charge  to  be  one,  and  r 
is  the  constant  electromagnetic  angular  momentum  of  the  particle. 
The  magnetic  field  B  may  be  determined  from  the  equation 


B  =  V  X  A 


(K2) 


with 


A  =  A(p,z)t. 


The  components  of  B  in  the  p  and  z  directions  are  then 


B  =— 


& 

"3z 


B  = 


(O) 


We  make  the  following  assumptions  on  the  magnetic  field  Bl 


(i)  The  field  strength  is  a  convex  function  along  a  segment 
S  of  a  magnetic  field  line 


58 


5 

\ 

(ii)  The  potential  A(p,z)  is  an  analytic  function  of  p 
and.  z  in  a  neighborhood  T  of  S. 

(iii)  The  magnetic  field  is  unequal  to  zero  on  S. 


Under  the  assumptions  (i-iii)  we  shall  prove  the  existence  of 
quasi-periodic  motions  in  T. 

(b)  The  main  step  in  our  proof  is  to  transform  the  Hamiltonian 
(4.1)  into  a  Hamiltonian  of  the  same  form  as  (3.6).  This  is 
accomplished  with  the  aid  of  the  following  Lemma. 

Lemma;  The  magnetic  lines  of  force  are  given  by  the  equations 

pA(p,  z)  =  c  =  const.  (4.4) 

Proof;  The  vector 


is  perpendicular  to  the  curve 


A(p,z)  -  |  =  0. 


(4.5) 


(4.6) 


From  (4.5)  we  may  write 


B  — 


3zP 


+  £  ^pA^z* 


(4.7) 
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Hence, 


which  vanishes  along  the  curve  A(P,  z)  =  c.  Thus,  r?*B  =  0,  and 
the  curves  (4.4)  coincide  with  the  magnetic  field  lines. 

Since  by  (iii)  the  magnetic  field  does  not  vanish  in  a  neighbor- 
> ood  of  S,  we  may  introduce  orthogonal  coordinates  a(P, z),  b(P, z) 
such  that  the  curves  a(p, z)  =  const,  define  a  magnetic  line  of 
force.  By  the  preceding  lemma  we  may  take 

a(p,z)  s  pA(p, z).  (4.8) 


The  coordinate  b(P, z)  is  then  determined  from  the  equation 

(4.9) 

To  effect  a  canonical  transformation  we  use  the  generating  function 


3a  3b  3h  3b  A 


F(a,b,pa,Pb)  =  a(p,z)pn  +  b(p,z)ph  , 


(4.10) 


and  the  Hamiltonian  (4.1)  assumes  the  form 


»(*,*,  Pa,  P„>  >  +  ^ 3  (h- “•) 


where 


,3ux2  <c&v2 

(■ss)  +  W 


,3bN2  ,c*k2 

<-3e>  +  <-3-> 


(it.n)' 


and  the  coordinates  p,  z  are  assumed  to  be  expressed  in  terns  of 

a  aid  b.  As  mentioned  in  the  previous  section,  the  quantity 
2  2  2 

li^P^  is  v„  .  To  prove  tills,  note  that 


^PPP  +  Pz"^  1/  3a*  3a^ 

- E - p(-  -Sp+ 

,  3a  3n  \ 

!p  +  ^  pz} 

r  /  3n.\  2 

r<^  p  1 


by  virtue  of  (4,3)  and  (4.8).  But 


3a  3n  ,  3n  3b  3n  3b* 

”3p  pz  ’  T  pp  =  ^  TE  '  S  ^’pb  * 


Hence 


3n  ch  3b  3b, 


6.1 


2 


p.  i.2^.2  2  2 


W '  ¥ 


src$>‘ft>‘«  <l>  <$>  ♦ 


♦  (#)2(|)2  *  <#>2(|n 


r,abN2  .  3bN2  2 

t  <■$>  +  <-s>  ifb 


-  V 


2  2 


b  * 


(e)  The  Hamiltonian  (4. 11)  is  now  very  similar  to  (3.6) 
and  we  proceed  in  the  same  manner  as  in  Section  3.  Our  first  step 
is  to  consider  only  those  particles  with  angular  momentum  r  for 
which  the  magnetic  line  of  force 


a(p, z)  *  pA(p,z)  =  T 

coincides  with  the  field  line  £.  (Of  course,  if  the  magnetic 
field  is  convex  along  every  field  line,  then  no  restrictions  are 
placed  on  T. )  As  in  the  previous  section  we  consider  only  those 
trajectories  which  stay  near  the  guiding  field  line  a(p, z)  =  r, 
and  which  have  small  energy.  Thus,  we  take 


a  -  T  =  ea 


P 


a 


% 


b  n 


4 


3 


j 


j 


h,  “  €IV 
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with  | a),  |pq!>  |Ppl  bounded  and  €  sufficiently  small  so  that 
2  2  2 

h  ,  h^,  and  p  are  analytic  functions  of  a  and  b.  The  new 
Hamiltonian  has  the  form 


rr/  n  \  1  -T  L,  2  2  ,  2  2\  ol 

^  H  =  ^(hapa  +  ly^)  +  _g  . 


(4.12) 


To  approximate  H  by  an  integrable  Hamiltonian  we  write 


hf  =  aL(b)  +  <?((a-l)) 


h^  =  b1(b)+  (a-1)) 


(^.13) 


ig  -  e  (b)  +  (a-1))* 

P 


The  Hamiltonian  (4,12)  may  then  be  written  in  the  form 


d  -  hq  +  nx  (4,i4) 

where 


H 


a 

T 


1/2  °1  2.  bl  2 

r(Pa  +  y*  Q  )  + 


(4.14)* 


and  ■  e1?1(a,€6,pa,Pp,6),  with  17^  analytic  in  all  its  vari¬ 
ables.  The  arguments  of  a^b^,  and  are  b  »  cp. 


* 
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Cur  next  step  is  to  introduce  new  variables  a* ,p^  in  place 

2  2 

of  a.p  so  that  H  is  a  function  of  (a*)  +  (p' )  alone. 

’  a.  c  '  v  a 

This  is  accomplished  by  taking  the  generating  function 


F(a,(3,P^)  =  K1ap^  +  pp£ 


where  =  c^^a^. 
form 


In  terns  of  the  new  variables  H  assumes  the 

o 


Ho“ 


, _  (a’)2+(p’)2  \  2 

i  +  4(p;r 


(^.15) 


Setting  a’  =  V*®  sin  0,  p'  =  cos  0,  we  note  that  to 

0( 

first  order  in  €,  R  is  constant  and 


0  - 


Yi  * 


The  quantity  is  the  magnetic  moment  M,  and  is  the  cy¬ 
clotron  frequency  B  evaluated  along  the  magnetic  field  line  i. 
To  prove  this  recall  that 


,dav2  .  21 


&1  "  ('3j)  +  V^) 


a  -  r 


1 

C1  •  "S 

P 


a  »  r  * 


Hence 


6U 


3  1  ,9a,2  1  ,Sa,2 

"l  °  Vi  = 

■  ^  +  <1  +  I'' 

-  B2. 


2.  2  2 

Since  e  D..p0  =  Vh  to  lowest  order,  it  now  follows  immediately 
Ip 

that 


R  = 


2  2 
€  V1 
2B 


to  lowest  order. 

Again  we  change  the  time  scale  so  that  the  frequency  be- 
comes  1J  our  trajectories  are  then  the  zero  energy  solutions  of 
the  Hamiltonian 


F  -  — (H  -  h) 
u>. 


(MO 


where  h  is  the  constant  value  of  the  Hamiltonian  M.  The  lowest 
order  term  F  of  F  is  now  given  by 


F  -  P.  - 
o 


(‘‘.13) 
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where  the  arguments  of  b^  and  are  ep'. 

To  determine  the  motion  in  the  3’,p^  plane  to  lowest  order, 
i.e.  to  determine  the  motion  along  the  guiding  field  line  I,  we 
consider  the  level  curves 


Mp, 


'f 

e' 


(^.19) 


in  the  3',p^  plane.  Under  the  assumption  that  s  B  is  a 
convex  function,  these  curves  will  be  closed  for 


B  , 
min 


< 


B 

c  M  nmx 


(**.») 


where  B^n  and  are  defined  in  Figure  15  below,  and  E 

is  the  energy  of  the  particle.  To  show  this  simply  write  the 
curves  l  in  the  form 


f;  O2  n  iLfh+aah ). 

Hence  the  notion  along  d  will  fee  periodic  for  c,h  satisfying 
(l.ao),  Ihe  legion  outside  is  usually  referred  to  as  the 

"less  cone"  of  the  particle,  since  one  cannot  hope  to  trap 
:  Articles  in  that  region. 


I 

I 

I 


Figure  15. 

Finally,  we  introduce  action  angle  variables  J,o  in  place 
of  (for  c  satisfying  (4.20)).  The  Hamiltonian  F  now 

may  be  written  in  the  form 

F  =  R  +  c(ej)  +  €j,e,$,e)  (4.21) 

where  ?1  is  analytic  in  all  its  variables.  This  is  exactly  the 
situation  we  encountered  for  the  dipole  field.  It  remains  to 
check  the  non- degeneracy  condition 


6? 


c"(r)  /  0 

with  r  =  cJ.  This  is  something  which  must  be  checked  for  each 
magnetic  field.  Moser’s  Theorem  will  then  guarantee  quasi- 
periodic  motions  in  the  region  c"(r)  4  0.  These  motions  densely 
cover  two  dimensional  invariant  tori  in  the  four  dimensional  phase 
space  of  the  particle.  Moreover,  any  trajectory  starting  between 
two  such  tori  can  never  escape.  Hence  all  particles  which  are 
adiabatically  trapped  are,  in  fact,  rigorously  trapped  for  all 
time.  Note  also  that  the  argument  in  the  preceding  section  con¬ 
cerning  the  existence  of  infinitely  many  periodic  solutions  is 
carried  over  exactly  to  the  more  general  case. 

(d)  To  conclude  this  section,  we  would  like  to  show  the 
impossibility  of  trapping  charged  particles  in  a  "planar"  magnetic 
field,  i.e.  a  field  whose  magnitude  and  direction  do  not  depend 
on  the  coordinate  z.  A  particle  moving  in  such  a  field  can  be 
described  by  the  Hamiltonian 

H(x>y>z>Px,Py,Pz)  =  |[PX  Py  +  (Pz  -  A(x,y))2].  (4.23) 

The  quantity  is  a  constant  of  the  motion  since  <&/bz  =  0. 

We  may  now  apply  our  method  to  get  quasi-periodic  motions  in  the 
x  -  y  plane.  The  motion  in  the  z- direct ion  is  then  found  by 
integrating  the  equation 


Z  =  Pz  -  A(x(t),y(t)) 


(4.&) 


However,  even  if  x(t)  and  y(t)  are  quasi-periodic  functions, 

we  cannot  expect  z(t)  to  be  quasi-periodic.  In  fact,  the  mean 

value  of  A(x(t),y(t))  need  not  be  p  ,  in  which  case  the  motion 

z 

in  the  z-direction  will  definitely  be  unbounded. 
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fells  report  is  concerned  with  rigorously  establishing  the  existence  of 
quasi-per iodic  motions  for  a  charged  particle  in  a  magnetic  field.  The 
following  results  are  obtained: 

(i)  The  existence  of  quasi-periodic  motions  lying  near  the  equatorial 
plane  for  a  charged  particle  moving  under  the  influence  of  the  earth's 
magnetic  field  (considered  here  as  a  magnetic  dipole).  This  proves  the 
stability  of  the  periodic  "equatorial  orbits"  for  small  values  of  energy. 

(it)  The  existence  of  quasi-per iodic  motions  which  need  not  lie  near 
the  equator,  exhibiting  thereby  the  trapping  of  particles  as  observed  in 
the  Von  Allen  radiation  belt.  Moreover,  these  motions  may  penetrate 
arbitrarily  close  to  the  dipole,  s  result  which  was  somewhat  unexpected. 

(Hi)  The  existence  of  quasi-periodlc  motions  in  a  rotacionully 
symmetric  magnetic  "minor  fLeld". 
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